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ABSTRACT: We study the dynamics of gravitationally collapsing massive shells in AdS 
spacetime, and show in detail how one can determine extremal surfaces traversing them. 
The results are used to solve the time evolution of the holographic entanglement entropy 
in a strongly coupled dual conformal gauge theory, which is is seen to exhibit a regime of 
linear growth independent of the shape of the boundary entangling region and the equation 
of state of the shell. Our exact results are finally compared to those of two commonly 
used approximation schemes, the Vaidya metric and the quasistatic limit, whose respective 
regions of validity we quantitatively determine. 
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1 Introduction 

By now, holography has established itself as one of the main tools used to gain insights 
into the out-of-equilibrium dynamics of strongly coupled field theories. Mapping the pro¬ 
cess of thermalization into black hole formation in asymptotically Anti-de Sitter (AdS) 
spacetime, gauge/gravity [1-3] methods have already solved several outstanding problems 
motivated by both heavy ion and condensed matter physics that have long eluded solutions 
using traditional field theory techniques (for reviews, see e.g. [4-7]). This can be largely 
attributed to the absence of competition: Perturbative methods typically fail already at 
moderate couplings, while time-dependent quantum phenomena are outside the realm of 
lattice Monte-Carlo simulations. 
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Important recent advances in applied holography include a fully dynamical description 
of shock wave collisions in strongly coupled J\T = 4 Super Yang-Mills (SYM) theory [8-13] 
as well as extensive work on the evolution of entropy-like quantities such as the holographic 
entanglement entropy (HEE) [14-23]. At the same time, technical leaps have been taken in 
the incorporation of inhomogeneities and anisotropies in thermalization dynamics [24-28], 
the development of a formalism to evaluate out-of-equilibrium Green’s functions [29-32], 
as well as the first studies of thermalization dynamics away from the infinite coupling limit 
[33-36] and in non-conformal backgrounds [37]. 

The above list of references clearly reflects an ongoing pursuit to take the holographic 
description of equilibration dynamics closer to the physical systems realized in nature, 
which are typically characterized by complicated initial states, finite coupling strength and 
N c , as well as broken conformal invariance. In this approach, one is typically confined to 
determining rather simple observables such as the temporal and spatial evolution of energy 
density or pressure. A different line of research concentrates on the simplest thermalization 
models available, but attempts to compute more complicated quantities, such as various 
off-equilibrium Green’s functions and other non-local observables. One prominent example 
of such models involves the gravitational collapse of an infinitesimally thin but massive 
shell in AdS space [38, 39]; following these papers, several works have addressed a variety of 
physical phenomena including particle production rates [40, 41], the chiral magnetic effect 
[42], jet quenching [43] and even elliptic flow [44], Most of these calculations, however, apply 
the so-called quasistatic approximation and assume the time scale related to the collapse 
to be parametrically larger than the other scales of interest, thus effectively considering the 
shell a static object [45]. 

In a preceding paper [46], we reported results from a set of calculations inspecting the 
falling shell model in a fully dynamical setup, where the shell follows a physical trajectory 
solved from the Einstein equations. The quantities considered in this context were the 
HEE and the Causal Holographic Information (CHI), which are both examples of geometric 
probes whose determination reduces to finding the area of some bulk hypersurface. As this 
involved rather complicated calculations requiring finding and matching extremal surfaces 
and geodesics in a time-dependent background, one of the aims of our current paper is to 
walk the reader through the technical details of this work. In addition, we will, however, 
present a considerably more thorough analysis of the HEE, comparing in particular its time 
evolution to results obtained in the quasistatic approximation and in the Vaidya metric. 
Here, we will find that during all times at least one of these approximation schemes is in 
a good quantitative agreement with the full results. We will also analyze the dynamics of 
the collapsing shell itself, and provide the full details of the construction of a coordinate 
system continuous at the location of the shell, briefly introduced already in [46]. 

In references [17, 18], it was noticed that in the Vaidya spacetime the entanglement 
entropy of large boundary regions exhibits linear increase in time for an extended period. 
The coefficient of this increase, ue, quantifies the rate, at which the time evolution entangles 
the subsystem to its surroundings. The authors of [17, 18] proposed an interesting conjecture 
that the value of ve computed for a collapse from AdS to the AdS-Schwarzschild spacetime 
might provide an upper bound for the rate of entanglement production in any relativistic 


- 2 - 


quantum field theory. Furthermore, it was argued here that the rate ue is a property of 
the final equilibrium state only, as it is only affected by the metric of the final black hole. 
One way of testing this proposal is to consider different initial states that evolve towards 
the same thermal state at late times — an exercise straightforwardly implementable in the 
collapsing shell model. As we will see, in all of our results the rate ue is indeed seen to be 
independent of the details of the shell trajectory, i.e. of the way the non-equilibrium initial 
state is prepared. Thus, we find evidence supporting the picture that ue is a property of 
the final equilibrium state only. 

Our paper is organized as follows: First, in section 2 and the corresponding appendices 
A, B and C, we provide technical details of our calculations, including solving for the 
shell dynamics, constructing a coordinate system that is continuous across the shell, and 
deriving continuity conditions for geodesics and extremal surfaces at the shell. In section 3, 
we then analyze the solutions to the shell equation of motion (EoM), while section 4 as well 
as appendix D are devoted to deriving the HEE and analyzing the corresponding results. 
Finally, in section 5 we compare our numerical findings to the quasistatic and Vaidya limits, 
analyzing the regions of validity of these approximation schemes, and in section 6 we draw 
our conclusions. 


2 Details of the calculation 


In this section, we introduce the machinery needed to obtain the time evolution of the HEE 
we are after. To this end, we first introduce our collapsing shell setup and derive the EoM 
of a shell falling in AdSs spacetime in section 2.1. Then, we derive a coordinate system 
continuous at the shell in section 2.2, which we use to write down junction conditions for 
extremal surfaces and more generic geometric probes intersecting the shell. Several details 
of the calculations are left to appendices A-C. 


2.1 Setup and shell dynamics 


Just as in [46], we work in a spacetime characterized by a negative cosmological constant, 
into which we immerse a thin massive shell, whose energy momentum tensor is proportional 
to a delta function in the radial coordinate. 1 Since both inside and outside the shell, the 
space is a solution to vacuum Einstein equations, we choose the inside metric to be that of 
an empty AdS Poincare patch and the outside metric the AdS-Schwarzschild solution with 
Schwarzschild radius rh, 


ds 2 


f±(r) 


-f±(r)dt 2 + 


dr 2 

/± 0 ) 


+ r 2 dx 2 



if r > r s 
if r <r s 


( 2 . 1 ) 

( 2 . 2 ) 


Here we have introduced a notation that we will be using throughout the calculation, 
where the subscripts + and — refer to quantities evalauted outside and inside the shell, 

1 Similarities of the thin shell setup and the fully back-reacted numerical solution of the Einstein-Klein- 

Gordon system are discussed in [47]. 
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respectively. It is important to note that although the metric functions f± themselves are 
time independent, the location of the shell r s . and thus the location of the discontinuity, 
are time dependent. The radial coordinate r and the spatial coordinates x are in addition 
assumed to be continuous at the shell. This means that there are two different and a priori 
unrelated time coordinates and t_, which we will later relate to each other. 

The coordinates on the shell are chosen to be the proper time of the shell and the 
spatial coordinates x, denoted by 

[f]=(r,x). (2.3) 

The embedding of the shell in the five-dimensional space is then given by 

[y M ] = (t s ±(r),r s (T),x) , (2.4) 


where p is an index running over the five coordinates of the AdSs space. Requiring that r 
is the proper time of the shell, we can further relate t s ± and r s to each other by writing 

ds 2 = —dr 2 = —f± i 2 ± dr 2 + ^ dr 2 . (2-5) 

J± 

Thus, the derivatives of t s ± and r s with respect to the proper time of the shell — denoted 
here by dots — are related by 

f± is± = \/f± + r 2 s . (2.6) 

In appendix A, we derive the EoM of the shell, given by eq. (A. 7). To evaluate its right- 
hand side, we need to specify the energy momentum content of the shell in the appropriate 
coordinate system. To this end, we employ the perfect fluid form, 

S'*- 7 = (p + p)u l v? +P7*- 7 , (2-7) 


where u l the four-velocity of the fluid and 7 * J the induced metric on the shell. This is in 
fact the most general possible energy momentum tensor when imposing translational and 
rotational symmetry in the x directions. Since the time coordinate in the £ coordinate 
system is the proper time of the shell, the coordinate system is in the rest frame of the 
fluid, and thus u = (1,0). The two independent non-zero components of equation (A. 7) 
then become 


^d^ 

r s dr 


~~ (V/- +r 2 s ~ \fh Tff) 

r s (■ Vf-+r 2 s ~ Vf+ + rf) 


-8vt55 p , 
-8ng 5 pr 2 , 


( 2 . 8 ) 

(2.9) 


from which we can derive a scaling law for the energy density, 


d_ 

dr s 


(Prl) 



( 2 . 10 ) 


Considering the simple equation of state (EoS) p = cp, we finally obtain 


p (x r s - 3 ( 1+c ) , 
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so that we can define a constant of motion M satisfying 


' s 

When the shell is pressureless (c = 0), the constant M is directly related to the conserved 
mass of the shell. 

After inserting eq. (2.12) to (2.8), we get as the final EoM of the shell 


■ 2 = M 2 _ /_ + /+ (/_ - /+)M +6c 

Ts 4 rf +6c 2 4 M 2 

where /_ and / + are evaluated at the shell, r = r s . It is noteworthy that the functional 
forms of f± are at this point still arbitrary, and that this equation is first order in time 
derivatives. The latter fact implies that solving it requires only one initial condition, e.g. the 
value of r s at some known time r, while the initial velocity is encoded in the constant M. 
This equation can be interpreted as the non-linear generalization of the conservation of 
kinetic and potential energy in Newtonian mechanics. 

If we now insert the explicit forms of f± from equation (2.2), we obtain from the above 


= ~T% + 


2 r 2 

' s 


+ 


M 2 

,4+6c 


8„,6c 


4 r 


+ 


rtr 


h's 

AM 2 


(2.14) 


Using as the initial conditions r s (r = 
value of M as 


M‘ 2 = rg +6c 


0) = ro, r s (r = 0) = 0, this allows us to solve the 


V7+M - Vf-( r o) 


(2.15) 


or using the explicit form of f±. 


M 2 


2 r 


6(l+c) 

0 




(2.16) 


Together with the equation of motion (2.14), this determines how the shell falls as a function 
of its proper time. If one on the other hand wants to EoM of the shell in terms of the 
coordinate time, or possibly relate the discontinuous time coordinates on the two sides of 
the shell to each other, one has to further use the relation 


dt- _ _ f + / /_ +r 2 

dt+ t s+ f- Y /+ + r 2 ’ 

which applies at the shell. 


(2.17) 


2.2 The junction conditions 

In order to eventually determine the time evolution of the entanglement entropy in the 
boundary field theory, we must be able to solve minimal surfaces in the spacetime containing 
a moving shell. In particular, we need to know how to join the minimal surfaces across the 
shell, i.e. how they refract at the shell. As we will review in section 4, the determination of a 
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minimal surface can be phrased as a variational problem, where one extremizes a functional 
of the generic form 

S = JraWW'WaUrh (2-18) 

where d a x^{a) = dx tl {a) / dcr a , with a denoting some set of coordinates on the minimal 
surface and x^(a) encoding the embedding of the surface in the spacetime. In this section, 
we will work out the refraction conditions following from extremizing a generic functional 
of the form (2.18). Thus, the results we obtain can be applied to any geometric probes in 
the spacetime, such as geodesics, string worldsheets and minimal area surfaces. 

Varying the action of eq. (2.18) leads to equations of motion for x M (cr), the Euler- 
Lagrange equations, that involve first derivatives of the metric. As we are dealing with a 
metric that is discontinuous, these equations will have delta function contributions from 
the derivatives. One way to derive junction conditions for x^(a) would be to integrate the 
EoMs across these singularities; in our case, this is, however, difficult to apply in practice, 
so we will use a different method. Namely, we will in the following explicitly construct a 
coordinate system, where the metric is continuous at the position of the shell. Working 
within it, the EoMs will have no delta function singularities, and therefore the solution 
x^(a) and all its first derivatives dx fJ ’(a)/dcr a will be continuous across the shell. Then, 
to obtain the junction conditions in the original coordinate system, we simply perform a 
coordinate transformation back to the original coordinates, where the discontinuities in the 
derivatives reappear from discontinuities in the coordinate transformation. 

To explicitly construct the coordinate system described above, we choose the timelike 
coordinate to be the proper time of the shell, r. Correspondingly, the required spatial 
coordinate is chosen to be the proper physical distance from the shell normal to it, which 
we denote by A and use to define our time slicing. Thus, our coordinate transformation has 
the form 

(t±,r,x)( t, A,x), (2.19) 

where a complication, however, arises from the fact that the normal vector of the shell is 
only defined at its location. This implies that we need to parallel transport this vector to 
cover the other parts of the spacetime. Intuitively, we start from the shell and then head 
out in the direction of the normal vector, parallel transporting it according to 

V n n = 0 . (2.20) 

This requirement is clearly nothing but the geodesic equation, meaning that our new spatial 
coordinate is simply the physical distance from the shell along a spacelike geodesic normal 
to the shell at its location. 

In order to determine the metric in this continuous coordinate system as well as to 
obtain the desired junction conditions, we need to know how the coordinates in the different 
coordinate systems are related to each other. Instead of obtaining explicit expressions for 
the coordinate transformation, it is, however, sufficient to merely calculate the values of the 
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partial derivatives 2 


(dt 


( dt 


\dr ’ \d\ 


( dr 






<9r 


<9A 


( 2 . 21 ) 


at the shell. This exercise is performed in appendix B. 

We will now proceed to compute the first total derivatives dx^/da a in the outside 
patch, transform them to the new coordinate system, and then transform them further to 
the inside patch. Using the chain rule, we can write the necessary derivatives in the form 


dt 


d 


( dt .i 


77^7 = 777 *+^)^^ = 7W 777 + ( ~ 7 \ ) 777 


da 

dr + 

da a 


da' 

da c 


-r+ (t(<t),A(<t)) = 


V dr 
dr. 


dr 


da° 

dr 


dr ) A da a 


dt .i 


+ 


d A 
dr+ 
OX 


dX 


da a 
dX 
da a ’ 


( 2 . 22 ) 

(2.23) 


where we will now drop the index a from a a to simplify our notation. From these expres¬ 
sions, we then solve 


dX 

da 

dr 

da 


dr+ \ dt.+ ( \ dr+ 


9r ) , 


da 


\ 9 - J ; 


da 


(dr+\ fdt+\ _ (dr+\ (dt+\ 

\^Jx\ dX Jr \ 9X )r\^ Jx 


( 9r + 

l dX 


dt+ 

da 


dt+ 

~5X 


dr+ 

da 


dr+ 

~d\ 


dt+ 

dr 


(£)» 


<9t+ 

dX 


(2.24) 


(2.25) 


which, when evaluated at the shell using the partial derivatives calculated in appendix B, 
gives further 


dX _ _ . dt+ \J f s + + r 2 dr + 
da s da f s+ da ’ 

— = /I i r? dt+ i ' s dr+ 
da s+ s da f s+ da 


(2.26) 

(2.27) 


Next, we use the chain rule to express dt/da and dr/da in the inside patch, 


dt- _ (dt- A dr (dt- A dX 

da \ dr ) x da + \ dX ) T da ' 

dr_ (dr_ A dr ( <9r_ A dX 

da \ dr ) A da \ dX ) T da ’ 

which, evaluated again at the shell, produces 


(2.28) 

(2.29) 


dt- 

da 

dr- 

da 


V fs- + r 2 s dr f s dX 

fs- da f s _ da ’ 

■ dr j —- —r^dX 

r ^ + Vfs-+ri—. 


’ da 


da 


(2.30) 

(2.31) 


“Here we have introduced the notation (§f) c familiar from thermodynamics to keep in mind which 
parameter is held constant as the other one is varied. To make our expressions somewhat more compact, 
we have also suppressed the arguments of our functions: When using the coordinates r and A, r and t are 
functions of both of these variables, whereas the time and position of the shell, t s and r s are functions of r 
only. Furthermore, / is a function of r and thus of both r and A while f s = f(r s (T)). 
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Finally, we insert eqs. (2.26) and (2.27) into the above equations to get the junction condi¬ 
tions 


da 
dr _ 


r=r 3 


da 


r=r a 


dt + 

da 

dt + 

da 


Ps-Ps+ - r 2 s dr+ 
/- da 

r s {P s+ -p s _)+^± 


-=r„ /-/+ 


(Ps+ - Ids-) 


— (Ps+Ps- - rf) 

~=r 0 J + 


(2.32) 

(2.33) 


where p s ± = \Jf s ± + r 2 . As a consistency check, we verify that in the limit where the shell 
vanishes, /_ —>•/+, both of these relations become identities. Also, in the limit where the 
velocity of the shell approaches the speed of light f s —> oo, the junction conditions reduce 
to the ones previously found in the Vaidya spacetime, cf. e.g. [18]. 

Interestingly, the above matching conditions are valid in a space with an arbitrary 
dimensionality, and one only needs to modify the metric functions / + and /_ in eq. (2.2). 
Also, in Appendix C we show how the conditions get modified, if one takes as the starting 
point of the calculation a more generic metric, where the dt 2 and dr 2 components are a 
priori not related to each other. 


3 Properties of the shell motion 


In this section, we perform a systematic study of the solutions of the shell EoM for different 
values of the EoS parameter c, defined through p = cp. For brevity, we will here denote 
t s ±{r) by simply t±. 


3.1 Simple example: c = —1/3 

Let us start by considering in detail the case of c = —1/3, which exhibits the same qual¬ 
itative features as the more general cases studied later, but is computationally somewhat 
simpler. For this value of c, the equations of motion namely reduce to 


■ 2 2 , 2 ( r 0 
r s = —r s + r 0 


U = 


r s\/rf ~r A h 


r 4 _ „4 ’ 


of which we can solve the first one by direct integration, producing 


1 f r ° drr 


1 


7 r 2 / 

'0 Jr s 


1 - ( i 

1 r 0 


= - arccos| — 
ro‘ 


4 2 


r s \- 


or equivalently 


r s (r) = royj cos(2r). 


(3.1) 

(3.2) 


(3.3) 


(3.4) 


From here, we see that for small and negative t the shell heads towards the boundary, while 
at r = ro or r = 0 it turns around and collapses. At the proper time r = 7r/4, the shell 
reaches the singularity at r = 0. 
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Figure 1 . Left: Shell trajectories for c = —1/3 and for ro = (4, 6,10, 25, 50) (bottom to top). The 
dashed curve represents here an ingoing null geodesic starting from the boundary at t+ = 0. Right: 
Shell trajectories for ro = 3 and 5 and for c = 0, 0.3, 0.33 and 1/3 (from left to right). The units 
in both figures are chosen such that rh = 1. 


Many of the interesting features of the shell trajectories become apparent only once 
the trajectory is expressed in terms of the time coordinate 1. 1 _. Solving eq. (3.2) leads to an 
expression for t + in terms of elliptic integrals, which is not particularly illuminating. We 
will thus rather take a step back and solve the EoM for dr s /dt +, obtained by taking the 
ratio of eqs. (3.1) and (3.2), 


dr s 

dt + 



(3.5) 


Solving for i+ from here, we obtain 



(3.6) 


where we have defined the integration variable u = r/ro. A numerical integration of eq. (3.6) 
is shown in figure 1 (left). It is clearly seen from here that all trajectories asymptotically 
approach the horizon at r = r'h = 1 with the same exponential rate as a null geodesic, but 
that the approach towards the null geodesic becomes faster when ro is increased. 

Both of the above features can be understood from the integral of eq. (3.6). Near the 
horizon, it is dominated by its lower limit, where we can approximate 


1 f du 

Ar h Jr s /r 0 U - r h /r 0 


1 

4 r h 


log (r s 


r h ) + -, 


(3.7) 


leading to the relation 


r h + Ce~ 4rht + 


(3.8) 


i.e. a null geodesic near the horizon. The exact same thing happens when ro is taken 
towards the boundary with r s /ro fixed to a small number: The integral is again dominated 
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by the lower limit of integration, and we can approximate 


4 


1 f du 1 

— / — 2 + ... =-h , 

r 0 Jr 3 /r o U 1 r s 


(3.9) 


which gives 



(3.10) 


identified as a null geodesic near the boundary. As ro —> oo, a boundary observer thus sees 
the shell apparoching a null geodesic, implying that the whole spacetime for r <C ro is well 
approximated by the Vaidya limit. 


3.2 Generic EoS: — 1 < c < 1/3 

For the range —1 < c < 1/3, the shell trajectories share the same qualitative features as 
the above example c = —1/3; in particular, the they are always seen to approach a null 
geodesic when either ?’o/ r /i —> oo or r s —> r^. To demonstrate this, we work at the level of 
the EoM and show that it approaches the equation of a null geodesic, 


dr 

dt + 


~f +( r ). 


in these limits. 

For a general c, the shell equations of motion are given by 


= —r + —-l 
s 2 r 2 


M 2 


4rf +6c 


+ 


r 8 r 6c 

r h r s 

4 M 2 ’ 


t± = 


Vf± + rj 

f± 


(3.11) 


(3.12) 

(3.13) 


where we will first consider the limit of the shell approaching the horizon, r s —>• r^. In this 
case, /+ approaches zero, so we can approximate (3.13) as 


4 


| 4 | 

4 ' 


(3.14) 


Using this, we obtain 


dr s _ rg_ 
dt + i + 




(3.15) 


which clearly implies that for all initial data with ro > r^ the shell approaches the speed of 
light, as it approaches the horizon. 

Moving next to the case of 7’o/r^ —> oo, we substitute the integration constant 


M 


—r 

2 


4„3c—1 
h r 0 


(3.16) 


into eq. (3.12). This leads to 


r 2 

' a 


Joc-2 


— T 

1 a 


2rl 


+ 


16 r? 


6c+4 


+ 


„6c 

s 


r 6c—2 > 
r 0 


(3.17) 
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which shows that as long as c < 1/3 and the shell location r s is kept fixed as ro/r^ —> oo, 
r s —> oo in this limit. Thus, in eq. (3.13) we can use r1 S> f+(r s ) and approximate 


4 ~ 



(3.18) 


which again leads to a null geodesic solution for r s . 

One difference between the two limits we have been considering above is that as the 
starting point of the shell is sent to infinity, the shell trajectory approaches a null geodesic 
when viewed either from the outside of the shell, using the time coordinate t+, or from the 
interior of the shell, using the time coordinate t_. This is not the case when ?’o is kept fixed 
and the shell approaches the horizon. In this case, only the trajectory as viewed from the 
exterior approches a null geodesic, while from the interior point of view it usually does not. 


3.3 Generic EoS: c > 1/3 

Finally, we take a look at values of c greater or equal to the conformal value c = 1/3. Again, 
the analysis of the trajectories near the horizon goes through unchanged, as the argument 
presented in the previous subsection is independent of the value of c. Thus, the shell is seen 
to approach the speed of light also for c > 1/3. 

The case of large r$ is, however, very different. Here, the AdS spacetime can be seen to 
provide a harmonic potential that pulls the shell towards the center, seen in the first term 
of the right hand side of eq. (3.12). For c > 1/3, the last term in this equation wins the 
pull of AdS, and the shell gets repelled from the center, accelerating towards the boundary. 
Physically, this means that the pressure of the shell wins over the gravitational attraction 
towards the center of AdS. 

In the special case of c = 1/3, the last term of eq. (3.12) also scales as r^, which leads 
to the possibility that the “forces” cancel at large r s . This enables the shell to approach the 
center of AdS with a very small acceleration even when it starts from near the boundary, 
as the gravitational pull and pressure almost cancel each other. Quantitatively, this can be 
seen by determining the acceleration of the shell, 


r s 


1 d .o 

2// dr r " ’ 


(3.19) 


and then expanding it near the turning point. Differentiating the expression in eq. (2.14) 
and expanding it yields 


r 


S 


—4r 0 + 2 (1 + 3c) r 0 



ro) ■ 


(3.20) 


From this expression, we see that the c = ^ case is special, as for sufficiently large ro the 
acceleration vanishes to first order and the motion of the shell can be arbitrarily slow. In 
a spacetime of arbitrary dimenson d + 1 (in our case d = 4), the special value of c reads 
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4 Entanglement entropy 


Next, we move on to consider the covariant holographic entanglement entropy (HEE) [48], 
which is obtained by extremizing the area functional 


A 


dr a 



dxv dx v 
da a da b 


(4.1) 


with the condition that the bulk hypersurface ends on a predefined surface A, which resides 
on a constant time slice on the boundary. In the dual CFT, the (geometric) entanglement 
entropy of the region V, whose boundary A is, is then conjectured to be given by 

See = 4 ^ ' (4 ' 2) 


The time evolution of this quantity has been extensively analyzed in various equilibration 
scenarios since the original work of [48]; in particular, for studies in the Vaidya spacetime, 
see [14, 15, 20]. 

In the current section, we will consider the evolution of the HEE in the collapsing 
shell model, using the physical shell trajectories obtained in the previous section. First, we 
study a simple example shape for the boundary region, a strip of width L, and then derive 
some more generic results for arbitrary shapes. To supplement this discussion, the relevant 
equations of motion for the case of a spherical boundary region are derived in some detail 
in appendix D. 

In all of the calculations presented in this section, we work in the Eddington-Finkelstein 
coordinates, where instead of using the time coordinate t we employ the lightcone coordinate 
v±. defined by 


dv± 


dt± + 


dr 

f ±( r ) ' 


(4.3) 


In addition, we will switch to the bulk radial coordinate z = 1/ r , so that our bulk metric 
will be given by 

ds 2 = —g [— h(z, v)dv 2 — 2 dvdz + dx 2 ] . (4.4) 

Here, we have further defined 


h(z , v) = 1 — 6(y — v s (z))z 4 , 


(4.5) 


where v s (z) is the trajectory of the shell parametrized as a function of z. In this entire 
section, we set the Schwarzschild radius to unity, i.e. rh = 1/Zh = 1. 


4.1 Strip boundary region 

The interior of a strip on the boundary is defined as the region of space with x l E 
{—L/2, L/2), x 2 E (—T 2 / 2 , L 2 / 2 ), and x 3 E (—L 3 / 2 , L 3 / 2 ), where L 2 and L 3 will be 
sent to infinity at the end. In this case, we can clearly assume that the bulk extremal sur¬ 
face is invariant under translations in the x 2 and x 3 directions. Thus, we can parametrize 
the extremal surface using the coordinates 2 = z{x) and v = v(x), where x = x 1 , while 
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the surface is spread homogenously along the x 2 and x 3 directions. For a more thorough 
explanation of the setup, we refer the reader to [15, 20]. 

With the above definitions, the area functional under consideration becomes 


A = L 2 L 3 / dx 


Vb 

z(x) 3 


= L 2 L 3 / dxC, B = 1 — h(z(x), v(x)) v'(x ) 2 — 2 v'(x)z'(x). (4.6) 


Due to the translational invariance of the system, there is a conserved quantity 

dC . dC . 1 

— z' + —v' — 


dzJ 


dv' 




(4.7) 


which is indeed constant along the entire extremal surface. Its value can be fixed by 
evaluating it at the point (z*, u*, x*) where the surface turns around, 3 , i.e. z ; (x*) = v'{x *) = 
0 which quickly leads to the result 


VB = r-* 

V z 


(4.8) 


The metric is clearly independent of v everywhere except at the position of the shell, 
implying that there is also a second constant of motion, 


a c 

dv' 


hv' + z' 

z z VB 


= ~E, 


(4.9) 


which takes different values on the two sides of the shell. Using eq. (4.8), we obtain 


h±v' + z' = E± , 


(4.10) 


where we have redefined the constant as E± = z*E±, and denoted h- = 1 and h + = 1 — z 4 . 
Solving this equation for v' and plugging the result into eq. (4.8) finally leads us to 


z ' 2 = E\ 


+ h± 


z*y 


-1 


= H ± (z) 


(4.11) 


If the boundary separation L is sufficiently small, the extremal surface never reaches 
the shell and always stays in the black hole region, implying that the entanglement entropy 
stays thermal at all times. The precise value of L, above which the surface crosses the shell, 
clearly depends both on the trajectory of the shell and on the boundary time. In the rest 
of this section, we will only consider the interesting case, where L is large enough so that 
the surface crosses the shell in the beginning of the time evolution. 

We start by studying the equations of motion in the pure AdS region inside the shell, 
where the extremal surfaces always have a turning point with z' = v' = 0. It is easy to see 
that FT vanishes there, which implies that everywhere inside the shell we have 


(4.12) 


Written in terms of the f_ coordinate, this means that t'__ = 0, i.e. that the surface lies in a 
constant time slice inside the shell. Integrating eq. (4.11) is now a straighforward excercise, 

^Here we assume that the extremal surface is reflection symmetric around x = 0. 
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and one finds a two parameter family of solutions parameterized by the turning point z * and 
the value of the time coordinate there, u( 2 *). The latter of these parameters can, however, 
be further traded for the point z c where the extremal surface crosses the shell, so that the 
interior surface is parameterized by the pair ( 2 *, z c ). In the following, we will need the value 
of the derivative at the interior shell position z'_ = z'(x c ), which is given by (cf. eq. (4.11)) 

= -^3 \/- z% . (4.13) 

z c 

Next, we continue the extremal surfaces across the shell using the junction conditions 
of eqs. (2.32) and (2.33), with a = x. In our current coordinate system, these read 

1 z 

z' + = — [a_|_a_ + z s (a- — a + ) — zf\ z!_ 4 —|(a_ — ot + )v'_ , (4-14) 

z c z c 

v' + = ^ 1 [a + a- - z s (a- - a+) - z%] v'_ , (4.15) 

Z C Il +\ Z c) 

where z± and v± are the corresponding derivatives evaluated on the outside and inside of 
the shell, a± = 1 Jh±z ■? + if, and z s is the proper velocity of the shell. Inside the shell, we 
can further use eq. (4.12) to write v'_ = —z!_. which reduces the junction conditions to 

7T = Z ( z s) = 4 (a+«- - if) , 

z c 

v' 1 9 

-T = V( z s) = —2 h (~ \ [“+«- - z s (a __ - a + ) - zf\ . (4.16) 

Z_ Z c tl-\-yZ c ) 


As a sidenote, we remark that we have here introduced a notation, where the junction 
conditions are considered functions of the derivative terms; this is done in anticipation of 
the following section, where we will consider the effects of the quasistatic approximation 
where these derivatives are altogether ignored. 

From eq. (4.16), we see that the quantity z'_, which depends on 2 * and z c . determines 
the values of the derivatives z' + and v' + outside the shell. Thus the integration constant E + 
gets determined by 2 * and z c using eq. (4.10), and we can therefore denote E+ = E+(z*, z c ). 
Now the boundary quantities can also be straightforwardly calculated using eqs. (4.10) and 
(4.11), and in particular the length of the boundary interval becomes 




max 


dz 

7^+M' 


(4.17) 


Here, we have denoted by z max the maximal value that the coordinate 2 obtains along our 
extremal surface within the outside region. If z' + < 0, then z max = z c as the surface climbs 
monotonically up towards the boundary, while if zL > 0, then z max is the point at which 
the surface turns around outside the shell, to be determined from the condition 


H+(z max ) = 0. (4.18) 

The time, at which the surface reaches the boundary, is on the other hand given by 


r z max si -y 

t = v s (z c )+ — 

Jz c "■+ 


E 




- 1 


L 


+ w- 


dz 

h. 


E 




+ 1 


(4.19) 
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while the area of the extremal surface reads 


A — 2L 2 L 3 z 3 


f 

J Zr 


dz 




f 

J Zr 


dz 




+ 


f 


dz 


i V / ^ 


(4.20) 


As we can see from here, all boundary quantities have now been given implicitly in terms 
of two parameters: the turning point z* and the crossing location z c . 


4.1.1 Early time behavior 

At early times, right after the shell is released from rest, the geometry is close to being static, 
and we can work in an expansion around a static shell. The relevant extremal surfaces then 
lie close to constant t surfaces, making it appropriate to use the (z, t) coordinate system. In 
particular, the shell trajectory near the turning point can be written in terms of a proper 
acceleration a = £ s (0) as 

z s (t) = zo + ^ar 2 + 0 (t 3 ), (4.21) 

where a can be determined from the equation of motion of the shell. 

Expanding now the junction conditions to first order in powers of z s = ar + 0(r 2 ) gives 


= y/h^z'_, t' + = 


1 - 


z c h h 


+ Z s z'_ 


(4.22) 


The boundary length and the area of the extremal surface are again given by eqs. (4.17) 
and (4.20), while the boundary time becomes 


t = t s (z c ) + 



dz 


E+ 

h+(z)^H + (z) 


(4.23) 


Here, t s {z c ) is once again the shell trajectory, now parametrized in terms of 2 and evaluated 
at the point where the extremal surface crosses the shell, z = z c . In the following we will for 
simplicity denote t s (z c ) = t c . Finally, the proper time r can at early times be approximated 
by the proper time measured by an observer at rest at z = z c , 


T 


yJh-\~ (z c )t c /z c . 


(4.24) 


At this point, an important observation is that t'< is proportional to z s ~ or, which 
is small at early times. Being proportional to t\. E + is therefore also small, and we can 
expand eqs. (4.17), (4.20) and (4.23) in powers of E + and 5z = z c — zq = ar 2 /2. Assuming 
L to be large, so that z * has to be sizable as well, we obtain for L 

\ = + °( z * 3 ) + Ofe-V). (4.25) 

We see from here that z * varies in time only very slowly, as the time dependence is supressed 
by an overall factor z*‘ i . Thus, z* is fixed in terms of L. 

To first order in E+, the boundary time becomes now 


t — t c + 


rz o 

/ dz 

Jo 


E a 


h+y/Ho® 


+ 0(t 3 ), 


(4.26) 
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*^EE(t) — ^ ee(O) 



Figure 2. The time dependence of the HEE (solid blue curves), obtained through a numerical 
integration of the extremal surface equations of motion, compared to the early time analytic solution 
of eq. (4.29) (dashed red curves). Here Zq = 0.5 and L = 8, while c = (—1, 0,0.1, 0.3,1/3) (left to 
right). 


where Hq(z) = H + (z) — E'+. For large z*, the integral can be easily evaluated, leading to 


t 


t c 


(1 - v 7 h+(zo )) 2 

- —s- a 

24 


(4.27) 


where we have used _ 

+ (z c ){ 1 - \A+(^c 


and replaced z c by zq, which is allowed to leading order in r. Similarly, we can expand the 
area in eq. (4.20) in powers of E + and 5z, which leads to the change in the area equaling 


)) 3 , 
—£* at c 


(4.28) 



AA = A(t) - A(t = 0) 


y/h+(zo)(l - y/h+izo)) 2 
T 2 E 3 -e- at c 


(1 - v^+Oo)) 2 

- —5 - a 

24 


Combining finally the above results produces 


at 2 


1 ( 1 ~ aA +("°)) 2 „ 

1 95 a 


1 V h +( z o) f 1 - y/h+ 

A^4 = - V 5 - 

2 4 


(4.29) 


where we denote the area of the boundary entangling surface as Aqa = 2 L 2 L 3 . This formula 
nicely demonstrates the relation between the entanglement growth and acceleration at early 
times. The early time behavior obtained here is compared to a full numerical integration of 
eqs. (4.17), (4.19) and (4.20) in fig. 2, which shows an impressive agreement up to relatively 
large time scales. 


16 - 






























V(z )/zt -V'(z )/z* 



Figure 3. The potential in the effective particle problem. Note that in this figure the boundary 
resides at z = 0, while the black hole singularity lives at z = oo. 


4.1.2 Linear scaling 

For large values of L, the quadratic early time behaviour of the entanglement entropy is 
followed by a long regime of linear increase, where ASee ~ At. In the case of Vaidya 
collapse, the existence of this region was demonstrated in [18], where it was seen to emerge 
from extremal surfaces inside the horizon of the black hole. In what follows, our analysis 
is closely related to the Vaidya case, and to this end we refer the interested reader to [18] 
for more details. Our goal will be to provide a simple and hopefully intuitive picture of 
where and why the linear region appears in our setup, highlighting the main differences 
that arise due to the slower motion of the shell. The precise details of the shell motion are 
unimportant for what follows, so we will only use the fact that z s (v) is a monotonically 
increasing function of v and that the shell does not move faster than the speed of light. A 
key assumption in our calculation is that the relevant extremal surfaces at late times are 
those that pass through the black hole horizon, which can indeed be shown to be true by 
numerically constructing the relevant surfaces. 

The equation of motion of extremal surfaces in the black hole region, cf. eq. (4.11), can 
be written in a suggestive form 

Z ' 2 + V(z)=£, v(z) = z* (J^ - , £ = E\, (4.30) 

where we have neglected two terms subleading at large z*. In the following, we will think of 
this equation as an EoM for a non-relativistic particle moving in the potential V(z), with 
x interpreted as a fictitious time coordinate. This potential is plotted in fig. 3. 

The rule for constructing the extremal surface is as follows. As we saw previously, just 
outside the shell the derivatives z’ and v' are determined by the two parameters z * and z c 
through the junction conditions of eq. (4.16). This setup is clearly equivalent to a classical 
mechanics problem with a particle starting from z = z c with some fixed initial velocity 
z 1 = vq and with the requirement of having to end up at the boundary. As can be seen 
from figure 3, the potential has a maximum at some z = z m , and for z > z m the force felt 
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by the particle tries to pull it towards the singularity. The location of this point is given by 


V\z m ) = 0 =* = 3 1 / 4 . 


(4.31) 


There are clearly two ways to avoid the fall of the particle into the black hole singularity. 
The first one is to have it start from z c > z m and give it enough negative initial velocity 
to get over the potential barrier at z m . This is, however, not allowed by the junction 
conditions, which determine the initial velocity through 

v 0 = ~^A(z c ), A(z c ) = y/z^h(z c ) + i 2 yjzl + i 2 - z 2 . (4.32) 

z c 

For general real values of z s , the function A(z c ) takes negative values in the region z c > z c . 
where 1 < z c < 2 1 / 4 and where the upper bound is approached when we approach the Vaidya 
spacetime as z s —> oo. This means that A(z^) is negative and vo positive for 2 C > z m , which 
implies that the particle will unavoidably fall into the singularity. 

A second way to reach the boundary is to start from z < z m and choose the initial 
velocity to be either positive but sufficiently small so that the particle will not get over the 
potential barrier, or alternatively even negative. For reasons that will become clear in a 
moment it is the first case that turns out to be the relevant one for us. Then, the maximum 
value for the energy of the particle is 

S ma x = V(z m ) = -^^-, (4.33) 

^rrt 


so that E .f should be bounded by y/£ m ax ■ In order to make the initial velocity vq given 
by eq. (4.32) small, we must obviously have the coefficient A(z c ) be very small. As A(z c ) 
changes sign at z c . which falls in the interval (1 ,z m ), we can do this by choosing z c to be 
sufficiently close to z c . 

Since we want the value of v, at which the particle reaches the boundary, to be large, it 
should spend a long fictitious time x in the bulk. After the above considerations, it should 
now be clear how to arrange this: We should choose the initial velocity to be such that the 
particle almost reaches z m and then turns around. To quantify this statement, we expand 
the potential around z = z m , obtaining 

z ' 2 - u 2 (z - Zm ) 2 = -uj 2 5£, cz 2 = ~^V"(z m ), 5£ = (£ max - £)/u 2 , (4.34) 

and then integrate this expression to get 


(2 I" dy _ 1 w / £ max - £ \ 

w Jvse \Jy 2 - 5£ w° g \ w 2 ) 


(4.35) 


Here, we have used the integration variable y = z m — z and neglected the contribution of 
the upper integration limit, as it is subleading in the limit of interest, 5£ —> 0. The factor 
of 2 in front of the integral is due to the symmetricity of the trajectory. 

The value of the boundary time coordinate of the extremal surface is determined by 


z yf£ - z’ 


(4.36) 
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As z' is very small for most of the time, we can treat the right hand side of this equation 
as a constant, producing 


t « — 


V&n 


h(z m ) 


Ax. 


The area of the extremal surface is then given by 

. 4 L 2 L 3 zl f dy 


2L2-L32; 


zfLw 


Ivss yV - 5S 4 


Ax , 


(4.37) 


(4.38) 


from which we obtain, using Ax = —t and £ ma x = — 

V t'max 


A zz 2 L 2 L 3 ^ h ’j 3 ™) t. (4.39) 

Inspecting the obtained result, we clearly observe that the area increases linearly in 
time. Following [18] in defining ve = \/~h(z m )/zf n as well as in denoting the thermal 
entropy density by Sth = 1/(4Gjv) and the area of the boundary entangling region by 
Aqa = 2 L 2 L 3 , we can now write 


See ~ AgASthVEt, (4.40) 

where the numerical value of ve is found to be 

y/2 

Ve = ¥A~ 0.620403. (4.41) 

This is exactly the result obtained for the Vaidya case in [18]. The above discussion should 
make it clear that the result is independent of the shell trajectory and on the precise form 
of the junction conditions. 


4.2 Linear scaling for general shapes 

In ref. [18], the authors argue that in order to isolate the linear regime in the time evolution 
of the HEE, one can take a limit in which turning point of the geodesic is sent to infinity 
before considering the time to be large. In this limit, the time evolution of the extremal 
surface is mostly due to the region inside the black hole horizon but outside the shell 
(i.e. within the AdS Schwarzschild metric), and the surface furthermore moves very little in 
the x direction in comparison with the part in the pure AdS region. To obtain the leading 
contribution to the area of the surface, it is therefore sufficient to approximate it as moving 
only in the v direction. This way, the area functional becomes 

A = ^/„^- fc ( |) 2 - 2 |. (4A2) 

where Aqa is again the area of the boundary theory entangling region and we have param¬ 
eterized the surface with v = v(z), denoting v' = dv/dz. 

Within the black hole region, there is again a conserved “Hamiltonian” 

d L hv' + 1 . A . 

— =- . = = -E , (4.43) 

dv' z 3 ^/—h(v ') 2 — 2 v' 
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from which we can easily solve 


v' = - ± 


E 


E 2 + h 


- 1 


(4.44) 


where the ± correspond to two branches of solutions. Upon integration, this finally gives 


t = v 


m - r 

J Z r 


dz 

h 


E 


E 2 + E 


*)-[ 


dz 

h 


E 


- 1 


E 2 + E 


(4.45) 


Here, the trajectory of the shell appears only via v s (z c ), which is obtained by inverting the 
function z s (v). 

Next, let us concentrate on the dominant contribution to the above z integral that 
originates from the region near z = z ma x■ As the first two terms of the Taylor expansion 
of y f.E 2 + ^ aroud this point vanish, the value of t diverges logarithmically, signaling that 
these specific values of E and z max correspond to a critical surface. This way, we are lead 
to the two conditions 


E 2 + HZraax) = 


d h(z max ) 


dz. 


= 0 , 


(4.46) 


max ^max 


of which the second one can be solved for the value of z max and the first one for E, giving 

(4.47) 


= Zm. = 3 1 / 4 , E = \ — 


h(z m ) 


•'max — 


Using these results, the area becomes 

dz 1 


A = Aqa 


f 

J Z r 


= + Aqa 
h JO 


f 

JO 


dz 


1 


E 2 + A Jo z~ J E 2 + h 


and matching the main logarithmic contributions of the two integrals produces 

A ~ Aqa ywg J ■ 

JA/Z> rn 

Finally, we approximate E by its value at the critical surface and thereby obtain 


A ~ AgAVEt , ve = \ l ~ 


h(z n 


(4.48) 


(4.49) 


(4.50) 


a result in full agreement with that of [18]. 


5 Comparison with common approximation schemes 

The collapsing shell model of gauge theory thermalization has been extensively used not 
only in the context of studying entropies, but also to evaluate various correlation functions, 
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Figure 4. Left: The HEE evaluated for a sphere of radius R = 2 for z o = 0.2 and c = 0, 0.33, 1/3 
(from left to right). Shown here are the full results (solid red lines) as well as the quasistatic (black 
dashed lines) approximation thereof. Right: A zoom-in to the c = 0 case, with the Vaidya limit 
(blue dotted line) added in the plot. 


corresponding e.g. to the electromagnetic current operator or the energy momentum tensor 
on the field theory side [33-35, 40, 41, 45]. While entanglement entropy calculations such 
as [14, 15, 17, 18, 20] are typically performed in the Vaidya limit of a lightlike shell, the 
Green’s functions are usually determined in the opposite ‘quasistatic’ approximation, in 
which the shell is taken to be a static object when formulating the junction conditions for 
the corresponding bulk fields. Physically, this approximation amounts to assuming that the 
time scale associated with the collapse of the shell is considerably larger than any other 
time scales relevant for the system, such as the inverse energy scale 1 /oj of the two-point 
function considered. 

The calculations we have performed in this paper for the HEE allow us to make an 
interesting comparison between our ‘exact’ results, derived for shells following realistic 
trajectories and employing the full junction conditions, and the corresponding quasistatic 
and Vaidya limits thereof. To obtain the former limit, we simply set z s = 0 in the junction 
conditions of eq. (4.16), reducing them to the simple forms 


= Z( 0 ) = 


a+ct- 


= V (0) = - 


a+a- 

Zc h +{Zc) 


(5.1) 


where z'_ is given by eq. (4.13). 4 At the same time, the Vaidya result is available by 
merely replacing the shell trajectory by an ingoing lighlike geodesic. Naively, we expect the 
quasistatic approximation to be valid only at the earliest times, i.e. near the turning point 
of the shell, while the Vaidya limit should be approached at late times. 

In figure 4, we display the result of the above comparison for the HEE evaluated for 
a sphere of radius R = 2. Somewhat to our surprise, we observe from the left figure that 
independent of the value of c, the static approximation appears to work rather well until 

4 The physical timelike trajectories are, however, used elsewhere in the calculation. 
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Figure 5. A plot of the ratios Z(z s )/Z( 0) (upper part) V(z s )/V( 0) (lower part) as functions of 
the boundary time. In both cases, the three curves correspond to the three cases displayed in fig. 4, 
i.e. Zq = 0.2 and c = 0, 0.33, 1/3 (from left to right). 


fairly late times and only breaks down when the shell is very close to the horizon; to aid the 
comparison, we have marked with dashed vertical lines the boundary times, at which those 
extremal surfaces that intersect the shell very close to the horizon, at z s = 0.99^, anchor to 
the boundary. At this point, the deviation of the quasistatic result from the full one is still 
only at the 10% level. In the right figure, we take a closer look at the c = 0 case and include 
in the figure also the Vaidya limit. We observe that the Vaidya result in turn gives a very 
good approximation of the full one already at relatively early times, and in particular that 
a combination of the quasistatic and Vaidya curves approximates the physical behavior of 
the entropy to a few percent level at all times. It is tempting to speculate, whether this 
observation would generalize to other, more complicated observables as well. 

To understand the observed behavior, it is instructive to inspect, how the full matching 
conditions of eq. (4.16) relate to the static approximation of eq. (5.1). The result of this 
comparison is displayed in fig. 5, where we plot Z(z s )/Z( 0) in the upper and V(z s )/V( 0) in 
the lower part of the figure, both given as functions of the boundary time. Shown are three 
curves corresponding to the three different values of the shell initial data already inspected 
in fig. 4: c = 0, 0.3 and 1/3, with zo set to 0.2 in each case. Comparing to fig. 4, we observe 
that in all three cases the deviation of the quasistatic HEE from the full result can to a good 
accuracy be attributed to the growth of the derivative terms in the junction conditions. In 
particular, the onset of the rapid growth of Z(z s )/Z(0) in fig. 5 coincides very accurately 
with the point of time, when the quasistatic entropy starts to visibly deviate from the full 
result in fig. 4. 

Having seen that the success of the quasistatic approximation can be traced back to 
the junction conditions for extremal surfaces, it is natural to ask, to what extent one can 
understand the reason for the slow turning on of their derivative terms. To this end, we 
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Figure 6. Comparison of the Vaidya (blue curve) and quasistatic (purple curve) limits of the 
HEE at d = 2, 4 and 6 , evaluated for a strip with L = 8 . 


now switch to the r, t coordinate system, in which the matching conditions can be shown 
to take the forms 

^= T (r s ) = ^-(P--P + ), (5.2) 

r -± = R(r s ) = ±-((} + P_-r 2 s ), (5.3) 

with the quasistatic limit corresponding to 


T( 0) = 0, R{ 0) 



(5.4) 


and the Vaidya one to 

T(oo) = “T uT’ R(x) = lj w^ ■ (5 ' 5) 

Expanding these functions around r = oo in the case of a d-dimensional field theory, we 
obtain in the quasistatic approximation 

T( 0) = 0 R(0 ) = 1 - ^ + 0 (r" M ), (5.6) 

while the Vaidya counterparts of these results become 

°°) = ^+ 2 ’ = ( 5 - 

From here, we see that for r> 1 the two matching conditions quickly approach each other, 
and that in particular the asymptotic behavior of R is independent of r s . This explains, 
why the quasistatic entropy approximates the full results, and even the Vaidya limit, so 
well, when the shell is released from close to the boundary. 

Finally, we remark that it is clear from the above results that the quasistatic approx¬ 
imation should successively improve, as the number of spatial dimensions in the system is 
increased — a direct consequence of the d-dependence of the blackening factor in the AdS 
Schwarzschild metric, h(z ) = 1 — z d . This effect is indeed clearly seen in the three plots 
displayed in fig. 6 , where we compare the quasistatic and Vaidya limits of the HEE at d = 2, 
4 and 6 . In the quasistatic results, the trajectory of the shell is taken to be lightlike and 
start from the boundary. 
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6 Conclusions 


Many of the existing studies of holographic equilibration have been performed within highly 
simplified models, which however have the virtue of allowing the determination of rather 
complicated physical quantities. One prominent example of this is the description of black 
hole formation via the gravitational collapse of a thin shell of matter in AdS spacetime, 
typically motivated as resulting from a rapid quench in the dual field theory system. Up to 
very recently, the determination of physical observables in this model has, however, been 
only possible with further simplifications, such as approximating the spacetime by its Vaidya 
limit, corresponding to a lightlike shell, or in the quasistatic approximation where the shell 
is taken to move arbitrarily slowly when formulating the so-called junction conditions. In 
our previous paper [46], we took the first steps towards overcoming these limitations in the 
case of the Holographic Entanglement Entropy (HEE). In particular, we considered there 
the time evolution of the HEE in the background of a shell following its physical trajectory, 
varying both the equation of state and turning point of the shell. Doing so, we were able to 
verify the earlier conjecture of [17] concerning the existence of a linear regime in the time 
evolution of the HEE that only depends on the properties of the final state of the system. 

In the paper at hand, we have continued work in the direction of [46]. In particular, 
we have studied the universality of the early and late time behaviors of the evolution of the 
HEE, and provided the first ever comparison of a fully dynamical thermalization calculation 
with its Vaidya and quasistatic limits. As expected, we observed that the quasistatic results 
provide a good approximation of the full ones at early times, while the Vaidya limit is 
approached at late times. Similarly, the quality of the quasistatic approximation improves 
when the shell EoS is chosen to minimize the rate of the collapse, while the Vaidya limit 
works better when the shell is allowed to accelerate more quickly. What we, however, found 
surprising was the fact that for many shell EoSs one can identify a brief period of overlap 
between the regions of validity of the quasistatic and Vaidya approximations, which seems to 
be at odds with the opposite nature of the two limits. To this end, we performed a detailed 
investigation of the reason of this behavior, tracing it back to the form of the junction 
conditions for extremal surfaces at the shell and to the slow onset of derivative terms in 
them. This observation suggests that the behavior may well be of somewhat universal 
nature, extending beyond the HEE, and raises hopes that the finding may eventually be 
used to simplify the holographic determination of many other dynamical quantities. 

Apart from improving the general understanding of thermalization dynamics, the most 
important outcome of our present work is the introduction of new technical tools to aid 
dynamical calculations within the collapsing shell model. In particular, the construction of 
a coordinate system continuous at the location of the shell and the derivation of junction 
conditions for extremal surfaces penetrating the shell are results that we hope will find 
applications in many forthcoming works. 
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A Equation of motion for the shell 


In this first appendix, we derive the equation of motion of a thin shell undergoing gravita¬ 
tional collapse in AdSs spacetime. To simplify our expressions, we suppress here the indices 
± indicating whether we are inside or outside of the shell. All identities involving / or t 
can be seen to hold both inside and outside of the shell. 

The unit normal vector of the shell is easily seen to read 


n r = Vf + rl and n* = y-, (A.l) 

where we have chosen the vector to point towards increasing r, and / is evaluated at the 
location of the shell, r = r s . Using this, we see that the only nonzero components of the 
shell’s extrinsic curvature, 


I<ij = n 


a 


( d' 2 y a r a <9/ dy^ \ 

\<)E,'<)£,' dg d%> ) 


are 


K-r 


r + 

2/f + 7j 


and. K xx - Kyy — K zz 


-r s Vf + r 2 s . 


The induced metric on the shell is on the other hand given by 


(A.2) 


(A.3) 


Ifij — diV 11 djy g^v , 


(A.4) 


so that 

Jrr — 1 and ^/xx — ^iyy — ^izz — ^s • 

Thus, the trace of the curvature tensor reads 



(A.5) 


(A.6) 


To derive the EoM of the shell, we use the Isreal junction condition [49]. It states that 
the difference of the extrinsic curvature between the inside and outside of the shell is related 
to the energy-momentum content of the object through 


[AJji Yj A] — 8tt^5 Sjj , 


(A.7) 










where the square brackets denote the difference between the inside and outside, 


[O] = Oinside - ^outside • (A.8) 

The LHS of eq. (A. 7) is zero for the non-diagonal terms, while for the diagonal terms the 
expression inside the square bracket is given by 

K tt ~ 7r tK = - — \Jf + r 2 , (A.9) 

r s 

K xx - 'yxxK = : (rWf + rf) ■ (A. 10) 


B Determining the metric in the continuous coordinate system 


In this appendix, we determine the values of the partial derivatives (fy) A , (^y) etc. needed 
to construct a coordinate system continuous at the shell, cf. section 2.2. 

The equations determining spacelike geodesics are given by 


2 / 






f dr dt d 2 t 

f dcr dcr + da 2 

= o, 

(B.l) 

f dx\ 2 d 2 r 


(B.2) 

f {di) + d7‘ 

= o, 

2 dr dx d?x 
r da da da 2 

= 0, 

(B.3) 


where ’ denotes a derivative with respect to r, and a is the affine parameter of the geodesic, 
identified as the proper length a = A. The trajectories we are interested in are defined at 
a constant x, so that the equations can be integrated to give 


(R4) 


The value of the integration constant A can be fixed by requiring that the geodesic points 
in the direction of the normal vector of the shell at its location, 


dx» 

~d\ 


= , 
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A=0 


where n^ has the form 


[n M ] = (f a /f, yjf + f 2 s , 0) 
Evaluating eq. (B.6) at the shell, we easily get 
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dt 
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^ = V f( r s) + A 2 = \Zf{r s ) + r 2 s , 
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from which we see that A = r s . 

Since eqs. (B.4) - (B.5) apply for spacelike geodesics, which by definition have constant 
r, these relations immediately produce the two partial derivatives 


Isa 


f{r) 


and 



(B.10) 


To get the derivatives (§y) A and (§^) x , we on the other hand need to differentiate the 
integrals of the equations of motion. Equation (B.5) can be integrated to yield 



dr 

\Zf( r )+f2(r) ’ 


which upon differentiation w.r.t. r leads to the expression 
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By changing the integration variable in this equation to A and requiring that the expression 
vanishes (A and r are by definition independent variables) then gives 
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dX 


f + r. 
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Finally, to get (f£) A , we solve t as an integral of A from eq. (B.4) and then differentiate it 
w.r.t. r to get 

/ f)+ \ r^ r r f 1 / fir \ 

(B.14) 
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which quickly leads us to 
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Now that we have all the necessary partial derivates at hand, we can see what our 
desired metric looks like. Setting dx = 0, we get 


ds 2 = -f dt 2 + -y 


= dr 2 

+ 2dXdr 


-f( a A V + 4f£ 


\drj x f\dTy x 
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+ dX 2 
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where the AA-component can be simplified using the fact that A was defined as a proper 
length, 
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The non-diagonal part of the metric is on the other hand given by 

*»—/(£') fin +!fin fll) 


dr) x \d A 


/ \dr 




/ / hs 


(B.17) 


implying that the metric is everywhere diagonal. Finally, the rr-component is given by the 
(by construction continuous) function 


/ dt\ l f dr 
9tt = ~ J + "7 ~rr~ 


J \drj x f \drj x 
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It is important to note that although the result for the rr component of the metric looks 
complicated and has dependence on the function /, on the shell it is equal to just —1, 
independent of the functional form of f(r). This reflects the fact that r was defined as the 
proper time of the shell. 

C Generalized junction conditions 

In this appendix, our goal is to generalize the junction conditions to the case where there 
are two unknown functions in the metric, 


ds 2 = — f (r) dt 2 - 1—— + r 2 dx 2 . (C.l) 

g{r) 

Following the above treatment, we again first construct the continuous coordinate system, 
and then proceed to derive the junction conditions. 

We will again take the location of the shell to be parameterized by (t s (r), r s (r)); this 
time the relation between r s and t s , however, reads 

e-Hf't 

while the normal vector of the shell is given by 


--,V9 + r 2 , 0,... • 


The geodesic equations are now seen to take the forms 


frt 

V +t=o ’ 


28 - 



which — taking into account that the affine parameter is the proper length of a space-like 
geodesic — can be integrated to 


fi = A, 

" 2 = 9 ( 1 + f )' 

Requiring finally that the geodesic is normal to the shell, 

dx^ 

~dX 

we obtain for the constant A 


A=o 


V 9 s 

At this point, we can again read off the necessary partial derivatives at 
obtaining 
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Inserting these into the relations (cf. eqs. (2.24)-(2.25)) 


dX 

da 

dr 

da 


f \ dt 4 . / \ dr + 

\~3V) X ~d'* ~ ~d^ 

(dr+\ fdt+\ _ (dr+\ (dt+\ 

\dr) x \dx) r \ax) T \dr ) x 

( dr+ \ dt+ _ / dt+ \ dr+ 

y d\ J da \ d\ ) t da 

9r+\ / dt+\ _ f dr+\ /<9i+ 


dX 


dr J X V dT Jx 


dX 


we get at the location of the shell 
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Using the chain rule, we next express dt/da and dr/da inside of the shell as 
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which, when evaluated at the shell, produce 
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Combining finally all the above results, we obtain as the generalized junction conditions 


dt- 

da 

dr- 

dcr 



( 0 + - P-) , 


(C.19) 

(C.20) 


where j3± are defined as before and the expression is evaluated at the location of the shell. 


D Spherical boundary region 


In this last appendix, we provide some details for the computation of the HEE in the case, 
where the boundary surface has the form of a sphere of radius R. In this case, the extremal 
surface is independent of the angular coordinates due to rotational symmetry, and we can 
parametrize it as z = z(p), v = v(p), with p being the radial coordinate on the field theory 
side (i.e. on the boundary). The area functional becomes then 

A = 47t J dp P ) 3 \JB = 47t J dp C , (D. 

where B is the same quantity as in the strip case except that the derivatives z' and v' are 
derivatives with respect to p. 

Due to the explicit appearance of p in the area functional, there are fewer conserved 
quantities this time, making the spherical case slightly more complicated than the strip one. 
There is, however, still a partial time translational invariance away from the shell, which 
gives rise to the conservation law 


dC p 2 (hv' + z') 
dv' z z VB 


(D.2) 


with E again taking different values on the two sides of the shell. In the interior region (as¬ 
suming the boundary radius to be large enough so that the extremal surface passes through 
the shell) there is a turning point at p = 0, where z' and v' vanish, which immediately tells 
us that E = 0 in the interior, or 

v' = -z'. (D.3) 

Applying this identity in the Euler-Lagrange equation for z(p), we obtain 


z(pz" + 2 z! z + 2 z') + 3p(l + z' 2 ) = 0 . 


(D.4) 


This equation has a one parameter family of solutions 

z(p) = sjzl - p 2 , (D.5) 


labeled by the turning point z*, identified as the most general regular solution with a turning 
point at z = z*. For a second order equation, one should specify two initial conditions, 
which in our case are chosen as z(0) = z* and 2/(0) = 0. 
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Outside the shell, we use eq. (D.2) to solve for v' and plug this into the Euler-Lagrange 
equation for z, giving 


2(p l zh + E+z 7 )z" + 2p 3 h(3pz' 2 + 2zz') + z(p 3 z' 2 (—ph! -\-Az') + E 2 fz & h') + Qf) i h 2 = 0. (D.6) 

This equation needs to be solved numerically. Noting that the interior surface satisfies 

v' = — z', the junction condition for the derivatives at the position of the shell is again 

given by eq. (4.16), with the derivatives now understood as derivatives with respect to p. 

This way the value of the constant E + is fixed to be 

) _ 

2 v' + z' + ’ 

which completes our exercise. 
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